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Abstract 

In spring 2006 a special meteorological situation occurred in the European Arctic 
region giving record high levels of air pollution. The synoptic situation resulted in 
extensive transport of pollution predominantly from agricultural fires in Eastern 
Europe into the Arctic region and record high air-pollution levels were measured at 
the Zeppelin observatory at Ny-Alesund (78° 54'N, 11° 53' E) in the period from 25 
April to 12 May, In the present study we investigate the optical properties of the 
aerosols from this extreme event and we estimate the radiative forcing of this episode, 


We examine the aerosol optical properties from the source region and into the 
European Arctic and explore the evolution of the episode and the changes in the 
optical properties. A number of sites in Eastern Europe, Northern Scandinavia and 
Svalbard are included in the study. In addition to AOD measurements, we explored 
lidar measurements from Minsk, ALOMAR (Arctic Lidar Observatory for Middle 
Atmosphere Research at Andenes) and Ny-Alesund. For the AERONET sites 
included (Minsk, Toravere, Homsund) we have further studied the evolution of the 
aerosol size. Importantly, at Svalbard it is consistency between the AERONET 
measurements and calculations of single scattering albedo based on aerosol chemical 
composition. We have found strong agreement between the satellite daily MODIS 
AOD and the ground-based AOD observations, This agreement is crucial for the 
radiative forcing calculations. We calculate a strong negative radiative forcing for the 
most polluted days employing the analysed ground based data, MODIS AOD and a 
multi-stream model for radiative transfer of solar radiation. 
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Introduction 


In the investigations of climate change, aerosols are of vital interest, as they have a 

direct impact on the radiative balance by scattering of solar radiation and absorption 
of solar and thermal radiation. The dominating process depends on the absorption and 
scattering characteristics of aerosols, which is a function of their composition, shape, 
and phase. The direct effect of aerosols is still connected with large uncertainty 
despite the huge scientific focus during the last decades (IPCC, 2007). Bellouin et al. 
(2005) have estimated the global aerosols direct radiative forcing in 2002 based on 
satellite measurements. They calculate a clear sky direct radiative forcing (DRF) of- 
0.8 W m' 2 with a standard deviation (a) of ± 0.1, indicating a significantly stronger 
DRF than present model estimates. The forcing is strongest over land though also 
connected with the highest uncertainty in these areas, This is explained by the high 
mean aerosol optical depth (AOD) over land at 550 nm of 0.13 ± 0.02. A recent joint 
model study of the DRF since pre industrial times has provided new estimates based 
on nine global models with comprehensive aerosol modules (AeroCom) (Schulz et ah, 
2006). The AeroCom study gave a global DRF annual estimate of - 0.2 W m" 2 with a 
o ss ± 0.2. These two studies clearly emphasize the large uncertainty still connected 
with the direct radiative effect of anthropogenic aerosols in . the atmosphere. 

DRF requires special attention in the polar region due to the surface conditions. Ice 
and snow give rise to high albedos, whereas water has a low albedo. As a 
consequence, the importance of light-absorbing aerosols is even larger in the Arctic 
than elsewhere, as atmospheric absorption is enhanced by the high surface albedo of 
snow and ice. Furthermore, the albedo of snow and ice may be reduced by deposition 
of black carbon (BC) containing aerosols (Hansen and Nazarenko, 2004). The total 
effect of anthropogenic aerosols in this region is largely uncertain and few studies are 
available (Quinn et al, 2007). 

At present, local and regional anthropogenic sources are almost absent in the Arctic 
region. Nevertheless Arctic haze, commonly present in springtime, is a well-known 
result of long-range transport from mid-latitude sources in Russia, Europe and North 
America (Stohl et al, 2006; Quinn et al., 2007; Law and Stohl, 2007). In combination 
with transport, favourable meteorological conditions with strong inversion in late 


winter and spring and little precipitation result in the high aerosol levels (Quinn et al., 
2007 and references therein), while the dominating light-absorbing aerosols 
transported into the region are black carbon containing aerosols from fossil fuel 
combustion. A recent study demonstrated that aerosols from biomass burning (BB) 
are a more important source of BC than previously thought (Stohl et al., 2006). 

To quantify the aerosol influence on the atmospheric radiative balance, knowledge of 
the AOD and single scattering albedo (SSA) is required (Hansen et al., 1997). SSA is 
a measure of an aerosol’s scattering efficiency. A SSA of 1 indicates a fully scattering 
aerosol with a strong cooling effect. Low SSA values indicate a heating effect. 
Haywood and Shine (1995) showed that the critical turnover value of the SSA is 
largely dependent on the surface reflectance. In the Arctic region, which has a high 
surface reflectance the turnover value would be high, implying that even moderate 
levels of absorbing aerosols might heat the atmosphere in this region. In fact, Pueschel 
and Kinne (1995) propose aerosols with SSA values as high as 0.98 may still have a 
heating effect on the atmosphere in this region. 

In spring 2006 a special meteorological situation occurred in the European Arctic 
region, which resulted in record high level of pollution (Stohl et al„ 2007; Treffeisen 
et al, 2007). The monthly mean temperatures for January, April and May were the 
highest ever recorded. The synoptic situation resulted in extensive transport of 
pollution into the region, and record high air-pollution levels were measured at the 
Zeppelin observatory in Ny-Alesund. A detailed description of the transport and 
pollution levels in Ny-Alesund is presented in Stohl et al. (2007). The record levels 
were observed at the end of April and in the beginning of May, with maximum values 
of AOD of 0.68 at 500 nm on the 2 May. This is the highest value measured since tire 
beginning of the measurements in 1991. A comparable mean Arctic haze level is 
0.133 at 532 nm, and the spring background level is 0,071 for the period 1991-1999 
(Herber et al. 2002). According to Stohl et al. (2007), the origin of the pollution event 
was agricultural fires expanding to forest fires in Eastern Europe. The maximum daily 
mean concentration of PM 0 .7 estimated from differential mobility particle sizer 
(DMPS) measurements was 29 pg m* 3 on the 3 May. Ozone (0 3 ) and carbon 
monoxide (CO) reached values of 83 ppb and in excess of 250 ppb, respectively. 
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which is the highest ever measured at the observatory. In fact, the previously 
maximum hourly 0 3 concentration ever recorded since measurements started in 1989 
was 61 ppb, clearly illustrating the extremity of tills event 

In the present study, we investigate the optical properties of the aerosols from the 
extreme smoke transport event in spring 2006. We examine the aerosol optical 
properties from the sources region and into the Arctic in order to explore the evolution 
of the plume and possible changes in optical properties. A number of sites located in 
Eastern Europe, Northern Scandinavia and Svalbard allow for a regional 
characterization of the aerosol optical properties by means of remote sensing 
instruments. Additionally, we use daily MODIS AOD products when available. The 
use of satellite data allows for a more complete picture of the episode and its 
evolution in time and space. Based on tire collected data sets of AOD from ground- 
based and MODIS instruments, single scattering albedos, height profiles of the 
aerosols, and in situ observations of the aerosol chemical composition in Ny-Alesund, 
we estimate the regional direct radiative forcing (DRF) of the aerosols from the 
extreme episode and compare it to the DRF for typical aerosol spring levels in the 
Arctic. 

2. Instrumentation and measurements 

The location of the sites included in this study is shown Figure 1. The information 
about the instruments and available data products is summarized in Table 1, We have 
collected data from three AERONET (Aerosol Robotic Network) (Holben et ah, 
1998) sites namely Homsund (Svalbard, Norway), Minsk (Belarus) and Toravere 
(Estonia). Spectral AOD and the AERONET inversion products from sun photometer 
measurements (size distribution, single scattering albedo, complex refractive index) 
are available for these sites (Dubovik et ah, 2000), The data are collected with the 
standard network Cimel sun photometer, which provides direct sun and sky radiance 
measurements at seven spectral channels in the range 340-1020 nm. The sites in 
Minsk and Toravere are both located close to the source region where the biomass 
burning occurred and in the pathway of the plume toward Scandinavia. Homsund is 
located on Spitsbergen (Svalbard), where the aerosol plume arrived after several days 
and more than 2000 km of atmospheric transport. 


Four more sites, Oslo, ALOMAR, Sodankylti and Ny-Alesund, providing ground- 
based AOD observations are located in the pathway of the plume. This allows for an 
investigation of the aerosol optical properties of the plume and how they are modified 
during the transport from the source region and into the Arctic. In Oslo and at 
ALOMAR Brewer spectroradiometers, designed for UV irradiance measurements and 
total ozone determination, are used to derive AOD at 320 nm (Cheymol et ah, 2006). 
We estimate the uncertainty in AOD to be ±0.05 based on comparison with a Cimel 
instrument at ALOMAR. This instrument is now permanently located at ALOMAR. 
Unfortunately, at the time period of the smoke episode the instrument was removed 
for calibration. Despite the limited accuracy of the Brewer AOD it allows a good 
tracking of the episode development at these two locations. 

At SodankylS Observatory, northern Finland (179 m asl), AOD measurements have 
been conducted by the Finnish Meteorological Institute (FMI) with a Precision Filter 
Radiometer (PFR) since 2004. AOD is determined at four spectral channels, as 
recommended by WMO; 368, 412, 500 and 862 nm (Wehrli, 2000; McArthur et ah, 
2003). SodankylS Observatory is part of the Pallas-Sodankylfi Global Atmosphere 
Watch (GAW) station. For details about the other activities at SodankylS, see G6mez- 
Amo et al. (2006). Complementary data used for interpretation of the SodankylS AOD 
retrievals are based on in situ measurements conducted at the Pallas GAW station, 
125 km NW of SodankylS. These include measurements of aerosol size number 
distribution and scattering coefficients. Hatakka et al. (2003), describe the 
measurements carried out at Pallas in detail. 

In Ny-Alesund the spectral AODs were measured with a sky-radiometer, Model 
POM-02 manufactured by Prede Co., Tokyo. The instrument was placed on the roof- 
top of the Rabben Station building in Ny-Alesund. The sky- radiometer includes 1 1 
channels at 315, 340, 380, 400, 500, 675, 870, 940, 1020, 1600 and 2200 nm. The 
AODs were calculated for 5 wavelengths at 400, 500, 675, 870 and 1020 nm. The 
accuracy of AODs obtained here is expected to be better than ± 0.005. 
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For Ny-Alesund we have also included the AOD measurements of the sun photometer 
SP1A at the AWIPEV station (The Alfred Wegener Institute for Polar and Marine 
Research) presented in Stohl et al. (2007). 

Clouds, especially thin and steady cirrus clouds, cause erroneously high AOD values 
hence cloud-affected measurements are screened from the data series. The AERONET 
(Smirnov et al., 2000) and PFR-GAW (Wehrli, 2000; Smirnov et al., 2000; Harrison 
et al., 1994) automatic cloud-screening algorithms have been used in this work. For 
the data from Ny-Alesund the cloud screening was based on separately sky radiance 
inversion analyses for aerosol optical properties. 

Three of the sites, Minsk, ALOMAR and Ny-Alesund, are also equipped with 
tropospheric lidars. The lidars in Minsk and at ALOMAR are both included in the 
European lidar network EARL1NET (Pappalardo et at., 2006). The Rayleigh-Mie- 
Raman lidar of the National Academy of Sciences of Belarus, Minsk, is also a part of 
the CIS-LiNet (lidar network over Former Soviet Union countries) since 2004 
(Chaikovsky et al., 2005). The Nd-YAG laser in Minsk emits radiation at 355, 532 
and 1064 nm. The five channel receiving system detects intensity of lidar signals at 
the wavelengths 355, 387 (N 2 Raman) 1064 nm, as well as co-polarized and cross- 
polarized components at the wavelength 532 nm. The vertical resolution of the signal 
is 15 m and the duration of normal measuring series is 200 s. 

ALOMAR is located at Andenes at the Andpya island outside the Atlantic coast of 
Norway, about 300 km north of the Arctic Circle. This tropospheric lidar started its 
operation in July 2005 and is the only sub-Arctic/Arctic site in the EARLINET 
(Pappalardo et al., 2006). The Rayleigh-Mie-Raman lidar is based on a Nd-YAG 
laser, which emits radiation at 355, 532 and 1064 nm. At present the system is 
equipped with five detector channels: 355, 387 nm (N 2 Raman), 532 nm cross- and 
perpendicular polarized, and 1064 nm. The nominal range and temporal resolution is 
7,5 m and 70 s, respectively. For more details see Frioud et al. (2006). 

A Micro-Pulse Lidar instalment was set up at Ny-Alesund in 1998 and upgraded 
through die MPLNBT project in 2002 (Welton et al., 2001). The National Institute of 
Polar Research, Japan operates the instrument at the AWIPEV stadon. The MPL 


employ a Nd/YLF laser, emitting radiation at wavelength of 523.5 nm. The data are 
stored at one-minute time resolutions and 30 m vertical resolution. Details regarding 
on-site maintenance, calibration techniques, description of the algorithm used and 
data products are given in Campbell et al. (2002), Shiobara et al. (2003, 2006). 

Measurements of aerosol chemical composition and their size distribution at Zeppelin 
observatory are used to calculate the single scattering albedo and die DRF. The 
Zeppelin observatory is situated at the Zeppelin Mountain (478 m asl) on the western 
coast of Svalbard (Norway). The Norwegian Institute for Polar Research maintains 
this observatory, while the Norwegian Institute for Air Research' (NILU) is 
responsible for the research acdvities at the site. The major inorganic anions (Cl', 
N0 3 \ SO4 2 ') and cadons (Ca 2+ , Mg 2+ K + , Na + , NH/) were extracted by water from 
aerosol filter samples collected on a daily basis using an open filter face NILU filter 
holder loaded with a 47 mm Teflon filter (Zefluor, 2 pm). The cations and anions 
were quantified using ion chromatography. NO3* and NH 4 + are both subject to positive 
and negative artefacts; here we use the sum of the particulate and the gas phase in our 
calculations. The ambient aerosol content of elemental carbon (EC) and organic 
carbon (OC) were obtained from aerosol filter samples, which were collected using a 
Leckel SEQ47/50, operating on a weekly basis. EC and OC were quantified by 
thermal-optical analysis, according to the NIOSH 5040 protocol (Birch and Cary, 
1996). The aerosol chemical properties are presented in Stolil et al. (2007). In the 
calculation of the single scattering albedo based on the chemical in situ measurements 
from Zeppelin we have used the chemical composition from Zeppelin as presented in 
Table 2. 

3. Aerosol optical properties In Northern Europe and the European 
Arctic during spring 2006 

3.1 Satellite aerosol optical depth during spring 2006 

To perform a regional analysis and explore the evolution of the aerosol plume in time 
and space, we have used Moderate Resolution Imaging Spectroradiometer (MODIS) 
on board Terra and Aqua. We have studied the region 50 - 85 °N, and 0-40 °E, 
which includes the source area and the transport pathway of the plume into the Arctic. 
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Figure 2 shows MODIS images at wavelengths in the visible region and how the 
smoke covers the northeast of Europe and the northern part of Scandinavia on two 
selected dates: 2 and 5 May. MODIS hot spot indicating the fires are represented by 
the red marks in the pictures. On 2 May, the smoke was transported across Finland, 
Sweden, traversing ALOMAR at the coast of Norway and up to Svalbard. On the 5 
May a plume went south of Norway and north along the Norwegian coast to Svalbard 
in the Arctic. 

In this section, we present the daily level 3 AOD at 550 nm retrieved from MODIS 
Collection Version 5 to give a regional overview of the AOD during the smoke 
transport event. We have combined the data from Aqua and Terra ((MYD08_M3.005 
and MOD08_M3.005). In the next section we present a more detailed analysis of the 
aerosol optical properties based on the ground-based observations made at the sites 
shown in Figure 1. 

In Polar Regions, the MODIS data coverage is low as satellite measurements of 
aerosol properties are difficult due to the conditions with high surface albedo, high 
cloud cover, large solar zenith angle, and thus a long path through the atmosphere. 
Further, the AOD retrievals are not performed under cloudy conditions or in regions 
with high surface albedo due to the snow and ice surface. To increase the data 
coverage we have combined the data from Terra and Aqua. The total data coverage in 
the region during the investigated period was 69%. 

MODIS daily AOD at 550 nm is shown in Figure 3 for the period 25 April - 12 May. 
The Figure shows high AOD in Eastern Europe from 25 April, and how the plume 
developed and expanded during the period, The MODIS AOD values in the centre of 
the plume were in the range from 0.8 to well above 1 for the entire period. During the 
first two days the centre of the plume was located in the Eastern Europe influencing 
the AERONET sites in Toravere and Minsk. Two days later AOD values ranging 
from 0.6-1 can be seen above the Barents Sea, and close to Svalbard. At this time the 
plume was already influencing the sites at Zeppelin, Ny-Alesund and Homsund. In 
the period 29 April - 1 May, the main aerosol plume covered Eastern Europe, southern 
parts of Finland, Sweden, and northern part of Norway influencing the ground-based 
sites in Minsk, Toravere, and ALOMAR. During 2-4 May, the centre of the plume 


was shifted more towards Finland, western part of Russia and northern parts of 
Sweden covering the site in SodankylH. During the period from 1 - 5 May, high AOD 
values were calculated in the Arctic region around Svalbard as well covering tire sites 
at Zeppelin and Homsund. On the 5 May particularly high AOD values were observed 
north of Svalbard with maximum values between 0.7-0.8. From 5-12 May AOD in 
the source region was gradually reduced and return to normal values from around 9 
May in Eastern Europe. In this period the plume was visible from satellite pictures 
along the coast of Norway, 

3.2 Optical properties of aerosols from ground based observations 
during spring 2006 

3.2.1 Spectral aerosol optical depth features 

The sun photometers located at the various sites included in this study (Figure 1) 
allow monitoring of the smoke episode under clear sky conditions. They provide 
higher temporal resolution and accuracy than the satellite retrievals, although with 
limited spatial information. 

As mentioned in the introduction, the levels on 2 and 3 May constitute the AOD 
record at Svalbard after more than 15 years of measurements. Table 3 provides a 
summary of the annual long-term mean AOD and the Angstrtim exponent for the 
Arctic sites, together with typical values during Arctic haze events, and the peak 
values during the spring 2006 episode (Herber et al., 2002; Toledano et al., 2006; 
Aaltonen et ai., 2006). The AOD summer background level is lower than the mean 
values indicated in this table (Tomasi et ai„ 2007). In general the AOD increases by 
30400 % during typical spring haze compared to the corresponding background 
levels. The smoke episode in May 2006 is exceptionally intense, with an AOD 3-4 
times higher than in typical Arctic haze periods. 

The time series of AOD at the Scandinavian sites is presented in Figure 4 a) and b). 
The AOD is shown in the left panel whereas the Angsfrtim exponent (a) is shown in 
the right panel. The episode over Scandinavia is clearly evident in the Figures. On 2 
May there is a peak in the AOD at all the Arctic sites: AODx^oo of 0.68 in Ny- 
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Alesund, 0.52 in Homsund, 0.B7 in SodankylS, and AOD **320 of 0.45 at ALOMAR. 
A first intense AOD peak was measured at ALOMAR on 30 April. The weather in 
SodankylU permits a more continuous monitoring of the episode. For this site we can 
see three peaks; the first elevated AOD peak on 27 April, which is also present in Ny- 
Alesund, and a second peak on 2 May. The third main peak observed on 3 May 
reached a maximum value of 0.87. 

The Brewer instrument maps the episode affecting Oslo. Unfortunately there are 
almost no data available from 22 April to 4 May due to cloudy conditions, so it is not 
possible to recognize the first plume taking place at the end of April. From 4 May, the 
AOD measurements indicated that the aerosol loading was high. After the peak on 7 
May (AODx= 32 o= 0.52) the AOD decreased slowly due to the stable conditions, and the 
background level was not reached until 14 May. 

The Angstrom exponents shown in Figure 4 were high, with values in the range from 
1.5-2 for the days when the episode was most intense. In Sodankyia, the daily mean 
a was 1.51 on 27 April, 1.72 on 2 May and 1.88 on 3 May, whereas more typical 
values were recorded after the event (1.01 on 10 May). In Homsund values of 1.77 on 
2 May and 1.92 on 3 May were recorded. The time evolution of a in Sodankyla shows 
an increase during the days when the episode was intense (1.53 on 27 April and 1.87 
on 3 May) compared to the less turbid days. The increase in a follows a similar 
pattern as the AOD, with two periods of enhanced levels centred on 27 of April and 2 
May (Figure 4 a)). This shows that the plume produced an increase in the slope of the 
spectral AOD, demonstrating that the fine mode fraction of the aerosols is increased. 
Tliis implies that the transported smoke aerosols were predominantly small, which is 
in agreement with the results of the SAFARI campaign and literature on biomass 
burning aerosols (Eck et al., 1999; O’Neill et ah, 2002; Reid et ah, 2005). The 
observed elevation of the AngstrOm exponent for aerosols from fires that are 
transported into the polar region is in agreement with the characterisation of Polar 
aerosols presented by Tomasi et al. (2007) and the results at Ny- Alesund for this 
particular episode reported by Treffeisen et al (2007). 


Data from the Arctic and sub-Arctic locations axe similar to the data obtained in the 
aerosol source region. This is shown in Figure 4 c) by means of the AOD and a 
measured at the ABRONET sites in Toravere and Minsk. (Note in this case tire 
different scales compared to Figure 4 a) and b).) The first remark is related to the 
intensity of the episode. The AOD levels in Estonia and Belarus were as high as 1.33 
in Toravere (29 April) and 1 .52 in Minsk (4 May). The long term annual mean AOD 
for the sites are 0.17 and 0.20, respectively, whereas the long-term mean Angstrom 
exponent is 1.30 and 1.28 (Table 3). However, for the period 15 April - 15 May 2006, 
the mean AOD (500 rnn) is 0.43 in Toravere and 0.44 in Minsk, For the 
corresponding a the means are 1.56 and 1.50, respectively. The episode is thus 
characterized by high AOD, with peaks almost one order of magnitude higher than the 
normal values, together with an increase in the Angstrdm exponent up to 1. 6-2.1 on 
days with enhanced aerosol loadings. 

The intensity of the episode changes with time, and a similar pattern for AOD in 
Toravere and Minsk is evident (except for 30 April). The two time series are in good 
agreement, indicating that aerosol plumes with similar characteristics are affecting 
both sites for almost one month. The agreement in the a parameter is also remarkable. 
The uncertainty in the Angstrbm exponent decreases when AOD increases. For the 
AOD measured in this period, the Angstrbm exponent is reliable, with less than 0.1 
absolute errors in a for the nominal AERONET absolute AOD errors, 0.01-0.02 
(Toledano et al, 2005). 

3.2.2 Size distributions and single scattering albedo 

The sky radiance measurements obtained with the Cimel sun photometer permit the 
retrieval of particle size distribution and complex refractive index (Dubovik and King, 
2000). Hence the optical parameters essential for the aerosol characterization, such as 
the single scattering albedo and the asymmetry parameter can be obtained. As 
indicated in Dubovik et al. (2002), the retrieval of the particle volume size distribution 
is adequate in practically all situations. However, the retrieved values of SSA and the 
complex refractive index must be utilized only from the retrievals obtained with high 
aerosol loading (AODx=440>0.5) and for solar zenith angle SZA>50°. For the three 
AERONET sites analysed, Homsund, Minsk and Toravere, the inversion products are 
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available. The radii, geometric standard deviations and volume fractions based on the 
inversions in die period 15 April - 15 May, is given in Table 4. The single scattering 
albedo at 440 nm and 1020 mil is also included based on inversions retrieved under 
the restricted conditions during the event. The number of restricted inversions is given 
in parenthesis. 

The volume size distributions for Minsk, Toravere and Homsund at selected dates and 
times are presented in the right panel of Figure 5a), b) and c), respectively. The left 
panels show the average volume size distribution based on all the available 
distributions from the smoke episode as a black curve. The red curve is a reference 
volume size distribution representing a spring situation with typical AOD values. 

The average median radii for the fine mode ill Toravere and Minsk are 0.16 pm, 
hence significantly larger than 0.12 pm in Hornsund. This is also the case with the 
coarse modes. However, the main difference between the size distributions for the 
sites close to the aerosol source region compared to Homsund, is the mode centred 
from 5-7 microns which is not present at Homsund. Compared with background size 
distributions, both the fine and the coarse mode concentrations are approximately 10 
times larger during the event in Minsk and Toravere. This is similar to the observed 
increase in the AOD during the smoke episode. Hence this indicates that both fine and 
coarse aerosols are emitted during the fires, or are produced by dynamic processes in 
the fresh smoke. However, at Homsund the fine mode concentration during the most 
intense days is 10 times higher than the reference day shown, whereas the coarse 
mode only increases by a factor of 1.5-2. The coarse aerosols observed close to the 
source are likely deposited during the Pansport to Svalbard, while smaller aerosols 
have longer atmospheric residence times and followed the plume into the Arctic 
region. 

In two cases (2 May and 3 May -17.00 UTC both days), with AOD x= 5 oo of 0.46 and 
0,43, the volume size distribution in Homsund (in Figure 5 c)) has values comparable 
to the size distributions in the source area with AOD j^joo around 0.9. According to 
Dubovik et al. (2002) errors up to 50% in the surface albedo do not have relevant 
impact on the aerosol size distribution. However, due to the special conditions in the 


Arctic region we have examined this further. We performed a set of simulations of the 
sky radiance using the UVSPEC radiative transfer model in order to estimate spectral 
diffuse down-welling global irradiance for different AOD, sun elevation and surface 
albedos. The simulations indicated that for AOD = 0.5, SZA « 60° and high 
surface albedo (snow) the diffuse irradiance increased by 45% at 440 nm compared to 
the same simulation for low surface albedo (water). For longer wavelengths, the 
increase is less (22% at 1020 nm). The effect is smaller as SZA increases (34% of 
increased diffuse irradiance at 44Gnm and 13% at 1020nm for SZA=75°). This could 
explain the differences in the size distributions retrieved in Homsund at large SZA 
(75°) compared to the observations around noon (SZA=61° on 2 May). 

Regarding the study of biomass burning aerosols in Dubovik et ai, (2002) their results 
indicate that typical single scattering albedo are ranging from 0.90 - 0.94 at 440 nm 
(Eck et al., 3999; Dubovik et ai., 2002), The values are expected to be somewhat 
lower at higher wavelengths. The SSA for Homsund, Minsk and Toravere from the 
AERONET inversions are shown in Figure 6 and the mean SSA is presented in Table 
4 for the period 15 April - 15 May. The mean SSA is 0.92 for Minsk and Toravere, 
and 0.98 for Homsund. In all cases SSA decreases with wavelength, with a mean of 
0.81, 0.88 and 0.96 at 1020 nm at the different locations, respectively. The SSA at the 
Zeppelin observatory in Ny-Alesund is calculated with Mie theory based on the 
chemical composition described in section 2 with the average size distribution as 
described above. Hygroscopic growth is taken into account in the Mie calculations 
and relative humidity data are taken from ECMWF reanalysis data for the actual time 
period. At the 2 May the SSA is calculated to 0.98 (range between 0.980-0.982) and 
for 3 May the corresponding values range between 0.974 and 0.981 at 440 nm. At 
1020 nm the calculated SSA is lower in accordance with the observations (modelled 
range between 0.880 and 0.928). Importantly, the calculated SSA is consistent with 
the retrieved SSA from the AERONET inversion in the Arctic region. 

It is apparent that the SSA increases with the distance from the source. This is in 
accordance with results from agricultural fires in southern Africa (Abel et al., 2003). 
This increase in SSA implies that the aging and transport of the aerosols from the 
plume and towards the Arctic region reduce the aerosol absorption, whereas the 
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scattering processes become more dominant H might be explained the deposition of 
the large aerosols during the transport. The large amount of coarse aerosols near the 
source (shown in Figure 5) most likely contains a fraction of absorbing aerosols. If 
these larger aerosols are deposited during the transport to the Arctic, the total aerosol 
absorption will be reduced. The chemical composition measured at Zeppelin resulted 
in high SSA and supports this explanation. Approximately 60% of the characterised 
aerosol chemical composition at Zeppelin was organic matter. This is remarkably high 
compared to a more typical situation (5-10%) (Stohl et al, 2007). Hence the relative 
contribution of the absorbing BC was small, even if the measured concentration was 
high. The common organic and inorganic components emitted from biomass burning 
are believed to have low absorption in the UV and visible region (Schkolnik et a!, 
2007; Myhre and Nielsen, 2004) thus scattering dominates the interaction with 
radiation leading to the high the SSA in this region. The influence and variations in 
tite relative humidity is also relevant. The relative humidity is much lower in Minsk 
and Toravere than in Homsund and Ny-Alesund. Sensitivity calculations show that 
employing the relative humidity observed in Minsk and Toravere, reduces die SSA at 
Ny-Alesund to values in the range of 0.95-0.96. The increase of SSA can also be a 
result of condensation and/or formation of secondary aerosols (Abel et al., 2003). 
Lack of chemical composition of the aerosols close to the source makes it difficult to 
quantify the different contributions to the observed changes in the optical properties 
and further explain the observations. A combination of different processes is the most 
plausible explanation. 

3,2.3 Comparison of ground based AOD and MODIS AOD observations 

The AOD from MOD1S aboard both the Terra and Aqua satellites is used in the 
radiative forcing calculations in section 4. Thus a comparison of these data and the 
available ground based AOD at the different locations is required. In the comparison 
we have used area average time series of the MODIS Atmosphere Daily Global 
Product. We have extracted ffxl 0 grid averaged MODIS data with the stations as the 
centre of the grids. SodankylS is one exception. Due to low data coverage during the 
event we decided to use 3°x3° grid averaged MODIS data when there are no data 
available above the site. We have compared the data to MODIS AOD extracted from 
l°xl° -grid when possible and the results were satisfactory. 


The results of the analysis of the time evolution of the episode are shown in Figure 7. 
The Figure shows the MODIS AQDx= 55 o from the daily overpasses of Aqua and Terra 
as the red and blue curves respectively in all panels. The satellite data are compared to 
hourly averaged ground-based data in the period 15 April - 12 May 2006. Despite the 
differences in the time resolution the results are in good agreement. In particular for 
Minsk and Toravere the agreement between the datasets is remarkably good. Yet, 
there are some peaks not captured by the MODIS data, which can be explained by the 
higher temporal resolution of the ground-based data. Further, the MODIS data is 
slightly lower than the ground based data for Minsk, This is expected as the ground- 
based data from Minsk is for 500 nm. 

Regarding SodankylS, Homsund and Ny-Alesund, the agreement is also good. Only 
the peak value observed 8 May at Homsund is not captured by the satellite 
measurements. However, this peak is not detected at Ny-Aiesund and it could be due 
to highcimis. 

AOD obtained from the Brewer instruments in Oslo and at ALOMAR are shown in 
the two mid panels. The AOD from the Brewers are for X=320 nm, which should 
result in higher AOD than for the MODIS data. Accordingly, the figures demonstrate 
that there are significant deviations in the absolute values of the AOD from die 
Brewers and from MODIS. AOD achieved from Brewer measurements is a relatively 
new method and it is connected with large uncertainty as it is highly influenced by the 
neutral density filter spectral transmittance (Cheymol et al. 2006). However, the 
purpose of the comparison between the Brewer AOD and the MODIS AOD is rather 
to manifest the time evolution of the AOD in Oslo and at Andenes. The results from 
ALOMAR and Oslo in Figure 7 dearly show that the elevated AOD levels detected 
by MODIS and obtained from the Brewers are analogous, and that the time series are 
parallel. Thus the transport of the aerosol plume into the Arctic region is clearly 
evident in both datasets and the AODs from the Brewer measurements justify the use 
of the satellite data in the calculation of the radiative forcing of the episode. 
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3.3 Profiles of optical properties of aerosols from ground based 
observations during spring 2006 

The development of the vertical distribution of the aerosol layer has been studied by 
tropospheric lidar measurements at three sites: Ny-Alesund, ALOMAR (Andenes), 
and Minsk. The lidar at the latter site is located dose to the source region. The 
corresponding profiles of the aerosol backscatter coefficients (ABC) are shown in 
Figure 8. 

For a detailed description about the transport and meteorological development during 
the episode see Stohl et al. (2006) and Treffeisen et at. (2007). Based on NOAA’s 
Hysplit model backtrajectories the air flow towards Minsk changes from southward 
(24 April) to more westward (27 April) and west/northwestward directions in 
beginning of May. On 24 April the aerosol plume from the fires filled the boundary 
layer in Minsk up to 1.8 km with an average backscatter coefficient of about 1.5-1 O' 6 
m 4 sr' 1 . On 27 April the aerosol load increases to 4*1 0‘ 6 m 4 sr 4 below 2 km. The 
maximum layer altitude increased to about 3 km in the beginning of May, with 
maximum ABCs of about 7-lff* m 4 sr* 1 and 4- 10"* m 4 sr' 1 on 3 and 4 May, 
respectively. Such high values are comparable to high aerosol loadings occurring 
above Minsk in the summer during easterly and southerly air flow, when air masses 
are brought towards the station from the Ukraine and Black Sea region (Wandinger et 
al, 2004). 

As observed by MODIS and the Brewer instrument at ALOMAR the plume reached 
the European Sub- Arctic region around 29 April. The ABC profiles from 25 April are 
characteristic for background aerosol load when air masses arrive at the site from 
sauth/westerly directions. On 1 May, ALOMAR is clearly influenced by the polluted 
air masses from the Eastern Europe with subsiding air masses coming from the 
south/east. The resulting ABC profile shows a maximum of about 1.6-10* m 4 sr* 1 
around 2 km altitude and a top layer altitude of -2.5 km. During the next few days the 
pollution reaching the station continue to increase, and the maximum of the ABC 
profile doubled at 3 May to -3.2 10"* m' 1 sr' 1 around 1.5 km, and the top layer altitude 
increased to 2.8 km. The highest ABC values at ALOMAR was 3.5-10'* m' 1 sr* 1 
around 2 km with a top layer altitude at -3.1 km. This was measured 5 May in the 


evening. On 6 May the air masses were coming from soulh/easterly direction bringing 
clear air, and the ABC profile measured by the lidar indicated clean background 
conditions. In the evening the same day again a higher aerosol load was observed, 
comparable to the layer seen on 1 May. These air masses were coming from Europe 
and were transported westward to the British Isle around the high-pressure system and 
northward. 

In Ny-Alesund the lidar profiles showed enhanced ABC values on 27 April. A 
maximum aerosol backscatter coefficient of 2.5-10'* m 4 sr' 1 at -1.5 km was measured 
on 2 May. Further, enhanced ABC values can be found up to about 4 km at that time 
at this Arctic site. However, the main aerosol layer is located below 2 km, which is 
about one kilometre lower than the top of the plume observed in Minsk. During the 
second part of 3 May and the first part of 4 May, a denser aerosol layer was observed 
aloft. While on 2 May back-trajectories between 1 and 4 km altitude show a 
downsloping trend, on 3 May only the 1 km altitude back- trajectory shows 
subsidence, while above 2 km a lifting of the air masses was observed. Due to the 
smoke episode the layer averaged aerosol-to-backscatter coefficient increased from 58 
±2 sr (20-21 April) to about 67 ±8 sr (2 May) and the spectral Angstrom coefficient 
increased from 1.4 ±0.1 to 1.7 ±0.1. This is consistent with the results presented in 
section 3.2,1; polluted air with a pronounced contribution of aerosols in the 
accumulation mode (see Figure 5). 

4. Radiative forcing of the May 2006 event 

Atmospheric aerosols are shown by observations alone or in combinations with model 
studies to have a global radiative impact. Further, there may be strong local effects in 
regions with high aerosol loadings (Haywood et al., 2003; Haywood et al., 1999; 
Kaufman et al, 2002; Myhre et al., 2007; Myhre et al, 2003a,b; Yu et al., 2006). The 
radiative effect of aerosols in the Arctic may differ from other regions due to the high 
surface albedo of snow and ice and during the summer period due to the high solar 
zenith angle that strengthen the radiative effect of scattering aerosols (Haywood and 
Shine, 1997; Myhre and Stordal, 2001). In the present study we have modelled the 
radiative effect of the aerosols using a multi-stream radiative transfer model (Myhre et 
al, 2007; 2003a) with die aerosols' chemical composition as described in Section 2, 
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the average volume size distribution described in Section 3.2.2, vertical profile of the 
aerosols based on the lidar observations in Section 3.3, and AOD from MODIS 
collection 5 data. Data for clouds and relative humidity are taken from ECMWF 
reanalysis for the actual time period. In our model calculations we calculate the AOD 
based on chemical composition and ambient relative humidity and thereafter this is 
scaled with the AOD from MODJS. In regions with missing AOD from MODIS we 
have used values from nearby grid points, 

Radiative forcing calculations are performed as a difference between the MODIS 
AOD for the smoke episode and a typical background aerosol layer with AOD of 0.05 
for tins time period of the year. Figure 9 shows regional cooling effects of the aerosols 
between 30 and 40 Wm' 2 in daily average for 2 May in the regions with highest AOD. 
It is also a large region with values between -10 and -30 Wm* 2 . For the 2 May the 
MODIS data show high AOD in the Barents Sea north of Russia as well and 
remarkable low forcings are calculated, around -35 Wm' 2 . The values for 3 May, 
which is shown in the lower panel, are slightly weaker and the strongest forcing has 
shifted to die east in accordance with the maximum AOD values. In comparison, a 
haze layer with AOD values typical for the season and the same optical properties will 
have maximum cooling effect around 5 Win -2 . Over Svalbard the radiative forcing is 
positive due to a much higher surface albedo than in the rest of the studied region. 
North of 80°N there are no MODIS data but sensitivity simulations show positive 
radiative forcing in this region although of smaller magnitude than shown in Figure 9. 
The forcings ate typically of less than 5 Wm* 2 but could reach 20 Wm* 2 in some few 
cloudy regions. Quinn et al, (2007) calculated a positive radiative forcing for typical 
haze in the Arctic region. They based their simulations on observed SSA, which was 
lower than what was observed under this smoke transport event. 

5. Conclusions 

We have investigated the optical properties and distribution of the aerosols produced 
during the agricultural fires in Eastern Europe in spring 2006, Based on the transport 
analysis of Stohl et al. (2007) and the available MODIS data for the period we have 
selected 6 different sites all influenced by the smoke plume during its transport into 
the Arctic region. We present the AOD foe all sites, and the Angstrdm exponent, the 


volume size distribution, single scattering albedo (SSA), and the vertical distribution 
of the aerosols when available. In addition, we have analysed the MODIS AOD for 
the region and time period and used all the collected data and information in regional 
radiative forcing calculations of this heavy smoke transport episode. 

The single scattering albedo in Minsk and Toravere is exclusively based on the 
AERONET inversion products. However, for the two sites at Svalbard it is based on 
both AERONET products from Homsund and in situ observations at Zeppelin, Ny- 
Alesund. Close to the source the SSA are low, 0.92 (at 440 nm) in both Minsk and 
Toravere. At Svalbard, the SSA is 0.98 (at 440 nm) for both sites, hence there are 
consistency between the AERONET products and the calculations based on the in situ 
aerosol chemical composition with regard to the single scattering albedo. 

The direct effect of aerosols in die Arctic region require special attention as die region 
is characterized by short wave forcing both during day and night in the summer, and 
high surface albedo. In the radiative forcing calculations we have used AOD from the 
MODIS instrument on board Aqua and Terra. The uncertainty in die MODIS AOD is 
larger than elsewhere as satellite measurements of aerosol properties in the Polar 
Regions are difficult due to high surface albedo, large solar zenith angle, and long 
path through the atmosphere. We have compared the satellite retrieved AOD to the 
available ground based AOD observations to verify the use of MODIS AOD in the 
radiative forcing calculations. The results are convincing, as we find high agreement 
at all sites. 

The radiative forcing of such an extreme smoke episode in the Arctic show dial the 
aerosols have a strong cooling effect above the ocean, and a much weaker heating 
effect above the ice and snow covered area. In total the scattering and thus the 
atmospheric cooling will be the dominating process when the SSA is as high as in this 
pollution event It has been proposed that increased haze events from more frequent 
forest fires in a warmer climate will enhance the Arctic warming (Law and Stohl, 
2007). However, based on our analysis this is not obvious. A central topic regarding 
the future is die surface albedo and its expected change. The climate effect of the 
aerosols in diis region is particularly sensitive to the surface albedo and the predicted 
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change in the ice and snow cover as well as the deposition of dark aerosols on ice and 
snow, will reduce the warming effect of the absorbing aerosols, 
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Tables 

Table l: Summary of sites, instruments and data products. 


Site 

Coordinates 

Instrument 

Products 

ALOMAR 

69° 16’ N 

Lidar 

AOD 320 nm, Backscatter profiles at 

(Andenes) 

16° 00’ E 

Brewer 

532 nm 

Homsund 

77° 00’ N 
15° 33’ E 

Cimel 

AOD (340-1020 nm), Angstrbm 
exponent, single scattering albedo 

Minsk 

53° 55’ N 
27° 36’ E 

Lidar 

Cimel 

AOD (4 channels 440- 1020 nm). 
Angstrom exponent, single scattering 
albedo, Backscatter profiles at 532 nm 




AOD (17 channels 350-1065 nm) 

Ny-Alesund 

78° 54’ N 
11° 53’ E 

Lidar 

SP1A 

Angstrom exponent, backscatterer 
profiles at 523.5 nm, chemical and 
physical properties 

Ny-Alesund 

78° 54’ N 
11° 53' E 

Prede 

AOD (1 1 channels 315-2200 nm), 
Angstrbm exponent 

Oslo 

59° 54 1 N 
10° 43* E 

Brewer 

AOD 320 nm 

SodankyiS 

67° 22’ N 
26° 37’ E 

PFR 

AOD (4 channels 368-862 nm), 
Angstrdm exponent 

Toravere 

58° 15’ N 
26° 27’ E 

Cimel 

AOD (4 channels 340-1020 nm), 
AngstrSm exponent, single scattering 
albedo 
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Table 2: Chemical composition at Zeppelin during 30 April - 7 May 2006. Table 3: AODx» 3 «>™, AngstrOm exponent, a, mean values over the specified period 


Component 

Concentration 

(pgia*) 

and data during die May 2006 smoke transport event. Typical values for Arctic haze 
are also included for the Arctic and sub- Arc tic sites. 

Elemental carbonaceous matter (EM) (EM=ECxl.l) 

0,26 





Typical Arctic 





Site 

Yearly mean 





Organic carbon (OC) 

3.53 





haze level 



Organic matter (OM=OCx 1.8) 

6.35 


Period 

AOD 

a 

AOD 

O 

AOD 

a 

S0 4 2 ' 

1.44 

ALOMAR 

2002-2006 

0.086 

1.54 

0.13 

1.65 

0.45* 




Homsund 

2003-2006 

0.04- 






N0 3 - 

0.90 

Ny-Alestmd* 

1992-2006 

0.07 

>1.5 

-0.12 

1.88** 

0.68 

1.87 

nh 4 + 

0.64 

Sodankyla 

2004-2006 

0.08 

1.54 

*** 

*** 

0.87 

1.93 


Toravere 

2002-2006 

0.17 

1.30 

Not relevant 

1.32 

2.12 

Minsk 

2002-2006 

0.20 

1.28 

1.52 

1.92 


♦at X=320 nm 
♦♦Herber et at. (2002) 

♦♦♦There are no typical Arctic haze values for Sodankyla as this was the first spring with AOD 
measurements at the site 
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Table 4: Summary of the aerosol size distributions and single scattering albedo for the 
inversions realized under AOD^onra 0.5 and SZA>50°: mean Median Radius (pm) and 
geometric standard deviation (<*), SSA and number of inversion data (N) available. 
The number of inversions with AOD440nm>0.5, SZA>50° and >21 symmetrical angles 
are in brackets. 


Site 





Parameter 

Minsk 

Toravere 

Hornsund 

Ny-Alesund 

Tvl 

0,160 

0.159 

0.122 

ND 

<Ti 

0.420 

0.326 

0.372 

ND 

Volume fraction, mode 1 

0.040 

0.060 

0.547 


fv2 

1.956 

1.484 

1.804 

ND 

cr2 

0.487 

0.344 

0.320 

ND 

Volume fraction, mode 2 

0.079 

0.084 

0.453 


r V 3 

7.193 

5.660 

- 


03 

0.425 

0.434 

- 

■ 

Volume fraction, mode 3 

0.881 

0.856 

- 

- 

SSA^Onm 

0.917 

0.924 

0.984 

0.98* 

SSA|020n/n 

0.811 

0.878 

0.962 

0.88-0.93* 

N data 

114 (63) 

189(41) 

44(7) 

- 


* Retrieved from in situ measurements at the Zeppelin observatory 


Figures 

Figure li The map shows the location of the measurement sites included in die study. 

Figure 2 : The plume is clearly seen in the MOD1S pictures of northeast Europe atid 
Northern Scandinavia, The red dots in Eastern Europe indicate fires. 

Figure 3: Aerosol optical depth at 550 nm for the period 25 April-12 May 2006 from 
the MODIS Collection Version 5. The data from Aqua and Terra is combined in the 
plots. 

Figure 4 : Left column: Aerosol optical depth for (a) Ny-Alesund, Homsund and, 
Sodankyia, (b) ALOMAR and Oslo (c) Minsk , Toravere. Right column: Angstrbm 
exponent for Homsund and, Sodankyift in the upper panel and for Toravere and Minsk 
in the lower panel. 

Figure 5: Volume aerosol size distributions during the smoke transport event: (a) 
Minsk, (b) Toravere, (c) Homsund. Right panel: Average distribution of the most 
intense period together with typical distribution. Right panel: Selected distributions 
during the period. AOD is given at 500 nm. 

Figure 6: Column single scattering albedo retrieved from Cimel photometer data at 
Homsund, Toravere and Minsk. 

Figure 7: Comparison of I°xl 0 grid averaged MODIS AOD and AOD from the 
ground based stations in the period 15 April - 15 May 2006. The MODIS AOD from 
Aqua and Terra is the red and blue curves respectively in all panels. 

Figure 8: Height profiles of backscatter coefficients of aerosols from ground based 
observations at the sites in Ny-Alesund, Andenes (ALOMAR), and Minsk during 
spring 2006. 

Figure 9: Radiative forcing of the episode relative to an aerosol situation in the spring 
with AOD=0.05. 
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Figures 



Figure 2: The plume is clearly seen in the MODJS pictures of northeast Europe and Northern 
Scandinavia. The red dots in Eastern Europe indicate fires. 




Figure 3: Aerosol optical depth at 550 tun for the period 25 Aptil-12 May 2006 from the MODIS 
Collection Version 5. The data from Aqua and Terra is combined in the plots. 






Figure 4: Left column; Aerosol optical depth for (a } Ny-Alesund, Homsund and, Sodankylti, (b) 
ALOMAR and Oslo (c) Minsk , Toravere . Right column: Angstrtim exponent for Homsund and, 
Sodankylti in the upper panel and for Toravere and Minsk in the lower panel. 
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Figure 5: Volume aerosol size distributions during the haze event: (a) Minsk , (b) Toravere , 

( c ) Homstmd. Right panel: Average distribution of the most intense period together with 
typical distribution. Right panel: Selected distributions during the periode. AOD is given at 
500 tun . 
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Figure 6: Column single scattering albedo retrieved from Cimel photometer data at Hornsund, 
Toravere and Minsk. 
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Figure 7: Comparison of l°xl° grid averaged MODIS AOD and AOD from the ground based 
stations in the period IS April - 15 May 2006. The MODIS AOD from Aqua and Terra is the red 
and blue cur\>es respectively in all panels. 
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Figure 8: Height profiles of backscaiter coefficients of aerosols from ground based 
observations at the sites in Ny-Alesund, Andenes (ALOMAR), and Minsk during spring 2006 








